arXiv:1506.00554vl [astro-ph.EP] 1 Jun 2015 


Mon. Not. R. Astron. Soc. OOP. n~[ fT8E 2014) Printed 2 June 2015 (MN style file v2.2) 


On the calibration of the relation between geometric 
albedo and polarimetric properties for the asteroids* 

A. Cellino^, S. Bagnulo^, R. Gil-Hutton^, P. Tanga^, M. Canada-Assandri^, 
and E. F. Tedesco^ 

^INAF - Osservatorio Astrofisico di Torino, 1-10025 Pino Torinese, Italy. E-mail: cellino@oato.inaf.it 
^Armagh Observatory, College Hill, Armagh BT61 9DG, UK. E-mail: sba@arm.ac.uk 

^ CASLEO and San Juan National University, San Juan, Argentina. E-mail: rgilhutton@casleo.gov.ar, micanada03@casleo.gov.ar 
Observatoire de la Cote d’Azur, Nice, France. E-mail: Paolo.Tanga@oca.eu 
^Planetary Science Institute, Tucson, USA. E-mail: eft@psi.edu 


Accepted 2015 May 22. Received 2015 May 20; in original form 2014 October 15 


ABSTRACT 

We present a new extensive analysis of the old problem of finding a satisfactory calibra¬ 
tion of the relation between the geometric albedo and some measurable polarization 
properties of the asteroids. To achieve our goals, we use all polarimetric data at our 
disposal. For the purposes of calibration, we use a limited sample of objects for which 
we can be conhdent to know the albedo with good accuracy, according to previous 
investigations of other authors. We find a new set of updated calibration coefficients 
for the classical slope - albedo relation, but we generalize our analysis and we consider 
also alternative possibilities, including the use of other polarimetric parameters, one 
being proposed here for the first time, and the possibility to exclude from best-ht 
analyzes the asteroids having low albedos. We also consider a possible parabolic fit of 


the whole set of data. 
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1 INTRODUCTION 

The geometric albedo of a planetary body illuminated by the 
Sun is the ratio of its brightness observed at zero phase an¬ 
gle (i.e., measured in conditions of ideal solar opposition 
to that of an idealized flat, Lambertiar0 disk having the 
same cross-section. The word “albedo” comes from the Latin 
word Albus, which means “white”. According to its defini¬ 
tion, therefore, the geometric albedo is the parameter used 
to indicate whether the surface of a given object illuminated 
by the Sun appears to be dark or bright. The albedo is wave¬ 
length dependent. In planetary science, the geometric albedo 

* Partly based on observations carried out at the Complejo As- 
tronomico El Leoncito, operated under agreement between the 
Consejo Nacional de Investigaciones Cientificas y Tecnicas de la 
Republica Argentina and the National Universities of La Plata, 
Cordoba, and San Juan. 

^ The phase angle is the angle between the directions to the Sun 
and to the observer as seen from the object 

^ That is, a surface having a luminous intensity directly propor¬ 
tional to the cosine of the angle between the observer’s line of sight 
and the surface normal (emission angle). A Lambertian surface 
exhibits a uniform radiance when viewed from any angle, because 
the projection of any given emitting area is also proportional to 
the cosine of the emission angle. 


minor planets, asteroids: general. 


has been traditionally measured in the standard Johnson V 
band centered around 0.55 /tm, and it is usually indicated 
in the literature using the symbol pv- 

The geometric albedo is a parameter of primary im¬ 
portance. Being an optical property of a sunlight-scattering 
surface, it must depend on composition, as well as on other 
properties characterizing the surface at different size scales, 
including macro- and micro-texture and porosity. All these 
properties are the product of the overall history of an ob¬ 
ject’s surface, and are determined by the interplay of phe¬ 
nomena as complex as collisions, local cratering, micro¬ 
seismology, space weathering, thermal phenomena, just to 
mention a few relevant processes. The fine structure and 
composition of the surface affects properties of the opti¬ 
cal emission which determine the results of many observ¬ 
ing techniques, including photometry and spectroscopy. It 
is particularly important in determining the state of po¬ 
larization of the scattered sunlight in different illumination 
conditions, this being the main subject of this paper. 

The geometric albedo is also a fundamental parameter 
when one wants to determine the size of a small Solar Sys¬ 
tem body, having at disposal photometric measurements at 
visible wavelengths. In particular, a measurement of bright¬ 
ness in V light is not sufficient to discriminate between a 
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large, dark object and a small, bright one, if the albedo is 
unknown. 

Geometric albedo should not be confused with the so- 
called Bond (or spherical) albedo. The Bond albedo is the 
fraction of incident sunlight that is scattered in all direc¬ 
tions and at all wavelengths. The Bond albedo is needed to 
estimate what fraction of incident radiation is actually ab¬ 
sorbed, and therefore contributes to the energy balance of 
the body determining its temperature. It is possible to define 
the Bond albedo at any given wavelenght A, e.g. , Av for that 
at the V-band of the Johnson UBV system, as A\ = q\ ■ p\ 
dMorrison and Lebofskvl Il979ll , where px is the geometric 
albedo at wavelength A, and qx is the so-called phase in¬ 
tegral, first defined bv iRussdl (Il916l l as the integral of the 
directionally scattered flux, integrated over all directions: 

pTT 

gA = 2 / <I>(A, q) sinada 

Jo 

where <l>(A,a) is the disk-integrated brightness of the ob¬ 
ject at phase angle a. Unfortunately, a determination of the 
phase integral, which requires in principle many measure¬ 
ments of the scattered sunlight at visible wavelengths ob¬ 
tained in different illumination conditions, is very hard to 
achieve, and is seldom available in practice. 

Asteroid sizes and albedos have been historically deter¬ 
mined mostly by means of measurements of the thermal flux 
at mid-IR wavelengths (the so-called thermal radiometry 
technique), generally using space-base d platforms like the 
IRAS and, more rece ntly, t he WISE dMasiero et al.l[201ll ~l 
and Akari dUsui et al.ll2013tl satellites. At thermal IR wave¬ 
lengths the received flux depends primarily on the size of the 
emitting objec10, and only weakly on the albedo. In particu¬ 
lar, the Bond albedo determines the fraction of the incident 
sunlight which is absorbed by the surface and is available 
to raise the temperature of the body. The temperature, in 
turn, determines the spectrum of the thermal emission in 
the IR. Since the Bond albedo of the asteroids is usually 
fairly low (in general well below 30%), most of the incident 
solar flux is actually absorbed by the body, and the inten¬ 
sity of the thermal flux turns out to be mostly dependent 
on the size, whereas the dependence upon relatively small 
differences in albedo is much weaker. Moreover, the com¬ 
putation of the geometric albedo from the Bond albedo, as 
mentioned above, would require a knowledge of the phase in¬ 
tegral, which is essentially unknown in the vast majority of 
cases. As a consequence, it is not really possible to solve si¬ 
multaneously for size and albedo in practical applications of 
the thermal radiometry technique. What is normally done is 
to derive the size from thermal IR data alone (assuming also 
that the objects have spherical shapes), and then determine 
the geometric albedo by using the known relation 

log(D) = 3.1236 - 0.2H - 0.5 log(pv) (1) 

where D is the diameter expressed in km (supposing that 
the object is spherical), H is the absolute magnitude and pv 
is the geometric albedo. To do so properly, at least some V 

® And on the temperature distribution across its surface, includ¬ 
ing also a contribution from the fraction of the body facing the 
observer but not illuminated by the Sun, when observing at non¬ 
zero phase angle 


magnitude measurements obtained during the same appari- 
tioqj in which the thermal flux of the object is measured, 
would be needed, in order to derive from them a reliable 
value of the absolute magnitude H. Unfortunately, in the 
real world no measurements of the V flux are really done 
in thermal radiometry campaigns, and H is directly taken 
from available catalogues. In turn, these H values are de¬ 
rived from V magnitude data (often of quite poor photo¬ 
metric quality), mostly obtained in different observing cir¬ 
cumstances, and using a photometric model of the variation 
of V magnitude as a function of pha se angle. As mentione d 
by several authors, (see, for instance. I Muinonen et al]l2010l l. 
the magnitude-phase relation for asteroids is described by 
means of parameters which are generally poorly known, and 
this introduces further errors in the geometric albedo deter¬ 
mination. 

In summary, it is difficult to obtain very accurate deter¬ 
minations of the geometric albedo of an asteroid based on 
thermal radiometry measurements. Based on the relation 
described by Eq. [T] the relative uncertainty on the albedo 
should be twice the relative error on the size, in ideal condi¬ 
tions. In practice, values as high as 50 % in geometric albedo, 
or even more for small and faint objects, are common even 
when the relative error on the size is of the order of 109ll3. 
The best results require measurements of the thermal flux to 
be obtained at different wavelengths in the thermal IR, an 
acceptable knowledge of the variation of V magnitude with 
phase, and detailed thermo-physical models which can be 
developed when a wealth of physical data is available from 
different observing techniques (including a knowledge of the 
shape and spin axis orientation). 

In this paper, we focus on another possible option to 
obtain estimates of the geometric albedo of asteroids, or 
other atmosphereless solar system bodies. This is based on 
measurements of the state of polarization of the sunlight 
scattered by the surface in different illumination conditions, 
and on the existence of empirical relations between geomet¬ 
ric albedo and polarimetric properties. Our present analysis 
is mainly devoted to summarize the state of the art of this 
application of asteroid polarimetry, and to provide one or 
more updated forms of the albedo - polarization relation¬ 
ship, sufficiently accurate to be used in practical applica¬ 
tions of asteroid polarimetry by the largest possible number 
of researchers in the future. We analyze the current obser¬ 
vational evidence taking into account an extensive data set 
available in the literatur^fl, including also observations car¬ 
ried out mostly at the Complejo Astronomico El Leoncito 
(San J uan, Argentina), that have been published only re¬ 
cently dGil Hutton et al]|2014l') . The derivation of the albedo 
from polarimetric data is a challenging problem which has 

The apparition of an asteroid is the interval of time (several 
weeks) before and after each solar opposition epoch, when the 
object becomes visible to the observers 

® This can be seen by plotting together for a comparison the 
albedos found for many thousands of asteroids observed by both 
the WISE and Akari satellites 

® Including the polarimetric data available at the NASA Plane¬ 
tary Data System at the URL address http://pds.jpl.nasa.gov/ 
(files maintained by D.F. Lupishko and I.N. Belskaya), and th e 
data published bviG il- Hu tton and Canada- Assandril 1120111 . 1^0121) : 
ICanada-Assandri al.rii2012l l 
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been open for a long time. In the present paper we take into 
account both traditional approaches as well as new possible 
developments suggested by the data at our disposal. 


2 ASTEROID POLARIMETRIC DATA 

Classical asteroid polarimetry consists of measurements of 
the linear polarization of the light received from asteroids 
observed at different phase angles. The observations give 
directly the degree of linear polarization and the position 
angle of the plane of polarization. This is usually measured 
with respect to the orientation of the direction perpendic¬ 
ular to the scattering plane, namely the plane containing 
the Sun, the observer and the target. According to elemen¬ 
tary physical considerations (Fresnel reflection) one should 
expect the scattered sunlight emerging from the surface of 
an atmosphereless planetary body to be linearly polarized 
along the direction perpendicular to the scattering plane. 
This expectation is only partly confirmed by the observa¬ 
tions. The asteroid light at visible wavelengths turns out to 
be, as expected, in a state of partial linear polarization, but 
in different observing circumstances the plane of linear po¬ 
larization is found to be either perpendicular (as expected) 
or parallel (and this is a priori unexpected) to the scatter¬ 
ing plane. It is therefore customary in asteroid polarimetry 
to express the degree of polarization as the ratio of the dif¬ 
ference of intensity of light beam component I± having the 
electric vector aligned along the plane perpendicular to the 
scattering plane minus the intensity 7[| of the component 
having the electric vector aligned parallel to that plane, di¬ 
vided by the sum of the two intensities. This parameter is 
usually indicated as Pi in the literature and is given by: 

P 

' ii± + h)’ 

According to its definition, the module of Pi is the de¬ 
gree of linear polarization of the received light as explained 
in elementary textbooks in physics (because 1± and /|| are 
found to be coincident with /max and 7min measured through 
a Polaroid), but the sign of Pi can be either positive or 
negative, depending on whether /max corresponds to I±, as 
should be expected based on elementary physics, a situation 
normally referred to as “positive polarization”. When /max 
is found to correspond to /y. Pi becomes negative, and this 
situation is called “negative polarization”. 

It is important to note that we will always use the Pi 
parameter throughout this paper every time we will refer 
to asteroid polarimetric measurements. Its value will always 
be expressed (and plotted in our figures) in percent (%). We 
also note that asteroids are not strongly polarized objects. 
The degree of linear polarization turns out to be usually 
below 2%. 

In asteroid polarimetry, by measuring Pi at differ¬ 
ent epochs, corresponding to different values of the phase 
angle, it is possible to obtain the so-called phase - po¬ 
larization curves. Some examples are shown in Fig. [T] 
Many data plotted in this Figure have been obtained dur¬ 
ing several observing campaigns carried out at CASLEO 
(Complejo AStronomico el LEOncito) in the province of 
San Juan (Argentina), using the 2.15 m Sahade telescope 
dGil Hutton et aDl2014l '). From Fig. [T1 it is easy to see that 



Phase angle (degrees) 


Phase angle (degrees) 


Figure 1. Examples of phase - polarization curves obtained for 
four large main belt asteroids belonging to different taxonomic 
classes: top left: the dwarf planet (1) Ceres (G-class); top right: 
(4) Vesta (V-class); Bottom left: (6) Hebe (S-class); bottom right: 
(21) Lutetia (for a long time considered to belong to the old M 
class, now included in the X complex). Open symbols refer to 
data available in the literature. Full symbols in the plots iden¬ 
tify observations obtained at the CASL EO observatory. Most o f 
them have been published only recently l|Gil Hutton et ^l2014^ . 
The best-fit lines correspond to the exponential-linear relation 
discussed in Sect. [5] 


asteroids belonging to very different taxonomic classes tend 
to exhibit phase - polarization curves which share in general 
terms a same kind of general morphology, but with differ¬ 
ences which can be easily seen, and re present some classica l 
results of asteroid polarimetry (see also lPenttila et ^l2005l ). 

In particular, all the curves are characterized by a “neg¬ 
ative polarization branch”, extending over an interval of 
phase angles between 0 degrees up to a value ao which 
is commonly found to be around 20° of phase (“inversion 
angle”). The extreme value of polarization in the negative 
branch is traditionally indicated as Pmin. Around the inver¬ 
sion angle the trend of Pi as a function of phase angle is 
mostly linear, and the slope of this linear increase is com¬ 
monly indicated as h. The interval of phase angles which is 
accessible to Earth-based observers extends in the best cases 
little over 30° when observing main belt asteroids (the possi¬ 
ble observing circumstances being determined by the orbital 
elements of the objects). The interval of possible phase an¬ 
gles extends up to much larger values in the case of objects 
which can be observed much closer to the Earth, as in the 
case of many near-Earth asteroids. 
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3 HOW TO CALIBRATE ANY POSSIBLE 

ALBEDO - POLARIZATION RELATION 

From Fig.[T] it is evident that, at the same phase angle, dif¬ 
ferent objects exhibit different degrees of polarization. Some 
objects are more polarized than others. This can be inter¬ 
preted in very general terms as bei ng a consequence of the 
classical Umov effect (|Umovl flflORl l. which states that the 
degree of polarization tends to be inversely proportional to 
the albedo, according also to laboratory experiments. 

The task of determining the albedo based on available 
phase-polarization curves is important in asteroid science. 
The foundatio ns have been laid long a go in some classi¬ 
cal pa pers, like lZellner and Gradiel (Il976ll and IZellner et all 
(Il977tl . Since then, several authors have tackled the same 
problem, with analyzes based on new sets of observations 
and/or laboratory experiments. The idea is to determine 
some suitable relation between the distinctive features of 
the phase - polarization curves of some selected sources, and 
their albedo, assumed to be known a priori with good accu¬ 
racy. The set of calibration objects must also be representa¬ 
tive of the whole population. 

We cannot assume a priori that all asteroids (which are 
a quite heterogeneous population), must respect a unique 
relation between albedo and polarization properties. Only a 
posteriori can we assess whether it is possible to accurately 
derive the albedo using a unique relationship independent 
on the taxonomic class. The main task is therefore to find 
a suitable representation of such a relation between albedo 
and polarization, to be tested over the maximum possible 
number of calibration objects belonging to different taxo¬ 
nomic classes. 

In this respect, one must analyze available polarimet- 
ric data for the largest possible sample of objects having a 
well-known albedo, in order to find evidence of some general 
and satisfactorily accurate relation between albedo and po- 
larimetric properties. Unfortunately, it is not so easy to im¬ 
plement this simple approach in practical terms. The use of 
laboratory experiments, based in particular on polarimetric 
measurements of meteorite samples, for calibration purposes 
seems to be natural, and was adopted in the 70s, but it is not 
exempt from problems. There are some technical difficulties, 
including the need of observing the specimens at zero phase 
angle to measure their albedos, somethi ng which is in g enera l 
not a trivial task. In a classical paper, IZellner et al.l ill977l l 
also mentioned some problems in correctly assessing the in¬ 
strumental polarization in the lab. Another general prob¬ 
lem is to ensure that the used meteorite samples are really 
representative of the behavior of asteroids as they appear 
to remote observers. Unfortunately, the meteorite samples 
have to be treated to reprod uce the texture of th eir asteroid 
parent bodies. For instance. I Zellner et al.l ill977l l noted that 
to find similarities between the phase-polarization curves of 
asteroids and meteorite samples, the latter have to be first 
crushed, and the way to do that has consequences on the de¬ 
rived polarization properties. Moreover, the authors noticed 
how difficult it is to eliminate from meteorite samples any 
source of terrestrial, post-impact alteration. For all these 
reasons, starting from the 90s most attempts of calibration 
of the albedo - polarization relation have been based on the 
direct use of asteroid albedo values obtained from other tech¬ 
niques of remote observation. 


Several papers have been based on the idea of using for 
calibration purposes some sets of asteroids for which the ge¬ 
ometric albedo had been derived from thermal radiometry 
observations, mostly consisting of old IRAS measurements, 
or, more recently, WISE data. This approach has the ad¬ 
vantage of being able to use for calibration many objects 
belonging to practically all known taxonomic classes. There 
are, however, several problems. First, albedos derived from 
thermal IR data are model-dependent. They depend on the 
choice of some parameters, which are needed to simulate 
the distribution of temperature on the asteroid surface, and 
the dependence of the irradiated thermal flux in different 
directions. Apart from a limited number of cases, albedo 
values determined by thermal radiometry data, particularly 
for objects for which we have little information coming from 
other sources, are simply too inaccurate for the purposes of 
a robust calibration. This is a consequence of the problems 
discussed in Section [T] concerning the general lack of simul¬ 
taneous photometric data at visible wavelengths, and the 
consequent use of values of absolute magnitude that are af¬ 
fected by large errors. To add some confusion, in the past the 
catalog of IRAS-based asteroid albedos changed with time. 
In p articular, the first publishe d catalog of IRAS albedo val¬ 
ues (iTedesco and Veederj|l99^ 1 had been built using thermal 
IR data coupled with estimated absolute magnitudes com¬ 
puted prior the introduction of t he (H, G) asteroid pho¬ 
tomet ric system. Using these data. iLupishko and MohamedI 
(Il996ll derived a first set of values for the calibration pa¬ 
rameters included in the so-called slope-albedo law (see Sec¬ 
tion mi. Subsequently, a new IRAS albedo catalog was pro¬ 
duced u sing absolute magn itudes expressed in the {H, G) 
system jTedesco et al.| l2002l l. and these values were used by 
ICellino et all ( 1999h to derive an alternative calibration of 
the slope-albedo law. This problem should be expected to 
arise again, due to the fact that lAU has recently recom- 
me nded the use of a new photometric system {H,Gi,G 2 , 
see iMuinonen et al.ll201oll . implying that the albedo cata¬ 
logues obtained using IRAS, and more recently, WISE data, 
should be updated again. 

Currently, different authors use different calibrations 
available in the literature, including, in addition to those 
just mention e d ab ove, also much ol d er ca librations by 
IZellner et al.l ill974l l and IZellner et al.l (Il977lj . which were 
mainly based on laboratory experiments. This is not an ideal 
situation, and one of the major goals of this paper is just 
to provide one or more updated forms of the albedo - polar¬ 
ization relationship, to be used by most researchers in the 
future, depending on the polarimetric data at their disposal. 

iMasiero et al.l (l2012l l proposed a new kind of calibration 
based on different polarimetric parameters and using a sam¬ 
ple of 177 asteroids for which the albedo has been estimated 
from WISE thermal radiometry data. In many cases, the 
objects were observed by WISE at fairly large phase angles, 
and the problem of assigning in these cases reliable values 
of corresponding magnitude in the visible is particularly dif¬ 
ficult. As a con sequence, in this pap er we will also make a 
new test of the IMasiero et akl (l2012| j approach, but using a 
different sample of objects having albedos not derived from 
ther mal radiometry data (see Section O. 

IShevchenko and Tedescol (1200611 proposed to use for cal¬ 
ibration purposes a limited sample of asteroids for which 
both the size is known with extremely high accuracy, and 
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also the absolute magnitude in the visible is well known, be¬ 
ing based on large data sets of available photometric data. 
As for the size, the most accurate values are certainly those 
obtained either in situ by space probes, or those obtained 
by accurate observations of stellar occultations. At least for 
some of these objects, also the absolute magnitude values 
listed in the catalogs can be reasonably reliable, although 
we should always remember that the absolute magnitude is 
not, strictly speaking, a fixed parameter, but it varies at 
different epochs, being dependent on the varying aspect an¬ 
gle of the object. The aspect angle determines the extent 
of the cross-section of the illuminated surface visible by the 
observer, and varies at different apparitions. This variation 
of visible cross-section depends on the overall shape of the 
object (it is zero for an ideal sphere) and on the orientation 
of the rotation axis. 

Limiting our analysis to the best-observed objects, for 
which the size and the absolute magnitude are supposed to 
be well known, the albedo pv can be derived using the rela¬ 
tion between size, albedo and absolute magnitude (Eq. [T]) 

The list of objects with r eliable albedo proposed by 
IShevchenko and Tedescol l|200(J l includes 61 objects. Among 
them, there are some of the largest and most observed aster¬ 
oids, including (1) Ceres, (2) Pallas, (3) Juno and (4) Vesta. 
We note that in the case of (4) Vesta, however, we use a 
different value of albedo, namely 0.35 ± 0.02, based on the 
most recent, and very accurate value of size measured in 
situ by the Dawn probe. The uncertainty in albedo for this 


asteroid depends on the fact that the disk-inte grated albedo 
tends to change at differe nt rotation angles (|Cellino_eLal] 
l2015l f . Many objects of the lShevchenko and Tedescol 1 20061 1 
list are much fainter (including some small targets of space 
missions) and several have never been observed in polarime- 
try. In this paper we f o llow the approach indicated by 
I Shevchenko and Tedescol ll2006f ). This does not mean that 
we are not aware of some probl ems: first, we know that the 
I Shevchenko and Tedescol (l2006li object list is now fairly old 
and needs an updating. This includes both considering a 
larger, currently available set of high-quality stellar occulta- 
tion data, as well as using more accurate values of the abso¬ 
lute magnitude, to be computed according to the new {H, 
Cl, G 2 ) photometric system adopted by the International 
Astronomical Union. We plan to produce an updated and 
possibly longer list of calibration targets in the near future, 
but we postpone this to a separate paper. In our present 
work we lay the foundations for any future analysis taking 
profit of a larger list of reliable asteroid albedos. The still 
limited data base of asteroid polarimetric observations is for 
the moment the main limiting factor for the investigations 
in this field. 

We have long been involved in an observing program 
of polarimetric observations o f asteroids belonging to the 
IShevchenko and Tedescol (l2006tl list, in order to improve sig¬ 
nificantly the coverage of the phase-polarization curves for 
these objects. So far, we have been able to obtain decent 
phase-polarization curves for only a limited sample of the 
whole list, taking also into account that some objects would 
require the availability of larger telescopes and/or better 
detectors, as well as a larger amount of dedicated observ¬ 
ing time. The results presented in this paper, which follow 
a pre vious preliminary analysis published by llCellino et al.l 
l2012|j . are already sufficient to find an updated set of cali- 
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Figure 2. The slop e-albedo relation, in log-log sca le, for 15 as¬ 
teroids of the list of IShevchenko and Tedescol ll20nfih for which a 
computation of the polarimetric slope is possible based on mea¬ 
surements obtained so far. Objects for which we have at disposal 
at least 10 polarimetric measurements are indicated by full, black 
symbols. Objects having a number of observations between 5 and 
10 are displayed using open, green symbols. The obtained polari¬ 
metric slopes of all the objects in this plot have been used in the 
computation of the linear best-fit that is plotted together with 
the individual data. The corresponding values of the Ci and C 2 
calibration coefficients are also indicated. 


bration parameters for the classical form of the slope-albedo 
law adopted by most authors in the past. The new data also 
allow us to explore new possible ways to express the relation 
between albedo and polarimetric properties, which will be 
probably the future in this field once the data-set of aster¬ 
oid polarimetric measurements will grow significantly in a 
hopefully not-too distant future. 


4 THE CLASSICAL APPROACH: THE 
SLOPE-ALBEDO “LAW” 

In practically all papers devoted in the past to this sub¬ 
ject, the relation between geometric albedo and polarimetric 
properties has been assumed to be one of the following ones: 


log(pv') = Cl log(/i) -1- C 2 

(2) 

log(pv) = C 3 log(Pmin) + C 4 

(3) 


In Eq. ([ 2 )), which was originally proposed as early as in the 
70s in the first pioneering investigations by B. Zellner and 
coworkers, h is the so-called polarimetric slope, namely the 
slope of the linear variation of Pi as a function of phase 
angle, measured at the inversion angle (see Section El). In 
Eq. ([3l), the polarimetric parameter is instead Pmin, namely 
the extreme value of negative polarization. 
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Most investiRations avai lable in the literature (see, for 
instance. ICellino et al.ll2012l . and references therein), have 
used Eq. which is normally known as the slope - albedo 
relation, or law. In fact, according to many authors, Eq. [3] 
leads to less (or much less) accurate results than Eq. (2] We 
will come back to this point in Section while in the rest 
of this Section we will focus on Eq. [2] 

The measurement of the polarimetric slope h should be 
done, in principle, by measuring the degree of linear polar¬ 
ization Pr in a narrow interval of phase angles surround¬ 
ing the inversion angle. In practical terms, however, the 
observers rarely have at disposal an ideal coverage of the 
phase - polarization curve, and often the polarimetric slope 
is derived by making a linear fit of a few Pr measurements, 
located not so close to the inversion angle as one would gen¬ 
erally hope. We will see in the next Sections some new pos¬ 
sible approach to derive h when one has at disposal a good 
coverage of the phase - polarization curve. For the moment, 
however, in a first tre atment of available polarimetric data 
for the objects of the IShevchenko and Tedescol (l2006fl list, 
we adopt the usual techniques, and we derive h from a lin¬ 
ear least-squares fit of all available Pr data. In particular, in 
order to use only homogeneous and high-quality data: 

• We limit our analysis to polarimetric measurements ob¬ 
tained in the standard V filter. 

• we only use values of linear polarization Pr having nom¬ 
inal errors less than 0.2%. 

• We use only polarimetric measurements obtained at 
phase angles larger than or equal to 14° of phase, a value 
generally well beyond the phase corresponding to Pmin, and 
in a region of the negative polarization branch where Pr 
starts to increase linearly with phase. 

• We require to have at least 5 accepted measurements, 
and that the interval of phase angles covered by the data is 
not less than 3°. 

A smaller number of measurements, and a narrower in¬ 
terval of covered phase angles, would make the adopted po¬ 
larimetric slopes more uncertain. The use of Pr data having 
fairly large error bars, up to 0.2, is suggested by the general 
scarcity of polarimetric data, but the effect of low-quality 
measurements is mitigated because in the computation of 
the best-fit curves, we weight the data according to the in¬ 
verse of the square of their associated errors . As f or the 
nominal errors of the I Shevchenko and Tedescol (l2006l ) albe¬ 
dos, which are not explicitly listed by the authors, we de¬ 
rived them using the quality codes listed in the above paper, 
according to their meaning as indicated by the authors. 

In this way, we were able to compute reliable polarimet¬ 
ric slopes for 15 cali bration objects. Taking the corr espond- 
ing albedos from the lShevchenko and Tedescol ll2006lj paper, 
we can plot h as a function of albedo, in a log-log scale, from 
which the coefficients Ci and C 2 in Eq. m can be derived 
using simple least-squares computations. The nominal errors 
of both the calibration albedos and of the derived polariza¬ 
tion slopes were taken into account in the computation. 

The results, including also the obtained values of the Ci 
and C 2 calibration parameters, are shown in Fig.[^ The val¬ 
ues of Cl and C 2 , together with their corresponding errors, 
are also given in Table 1. A few considerations are suggested 
by looking at Fig.[2] First, the linear best-fit solution seems, 
as a first approximation, fairly reasonable. This seems to 



Py (Shevchenko cSc Tedesco 2006) 


Figure 3. Comparison between the albedos of 15 objects ob¬ 
served in our campaign, derived using our new calibration coeffi¬ 
cients of the slo pe-albedo relation, and the corr esponding albedo 
values given bv fshevchenko and Tedescol {200^. The meaning of 
the symbols is the same as in Fig. [2] 
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Py (Shevchenko & Tedesco 2006) 


Figure 4. Plot of the difference between the albedo values derived 
using the new calibration co efficients of the slope-albedo rela tion 
for objects included in the IShevchenko and Tedesc'3 l(200^ list 
and the corresponding values found by the above authors. The 
dashed lines correspond to differences of ±0.04 in albedo. The 
meaning of the symbols is the same as in Fig. [2] 
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Pofarimetric Slope (%/deg) 

Figure 5. A comparison between the calibration of the slope- 
albedo relation pre sented in this paper (red line) and those by 
ICellino et alj ||2C)12| 1 (black, dashed line). 

be confirmed by Figure [21 where the albedo values of the 
objects considered in our analysis, as they can be derived 
from our updated determination of the calibration coeffi¬ 
cients, are plo tted versus the corresponding albedo values 
determined bv IShevchenko and Tedescol (l2006l) . In Fig.|4]we 
show for each object the difference between the albedo value 
obtained from the polarimetric slope and the albedo value 
given by [Shevchenko and Tedescol (l2006fl . From Figs. [2] to 
[4] it can be seen that the discrepancies are generally low 
in absolute terms, being mostly below ±0.04, as shown in 
Fig-Hl The error bars of the obtained albedo values tend to 
increase with albedo, but this should be expected, because 
the slope-albedo relation (Eq.[2J implies that the error of pv 
must increase linearly with pv itselflll 

One object turns out to have a polarimetrically- 
derived albedo that is significantly discrepant, namely (2) 
Pallas, (the point which is located at the highest ver- 
tical distance above the bes t-fit line in Fig. ID). The 
[Shevchenko and Tedescol (l2006l ) albedo value of this aster¬ 
oid, 0.145, seems to be noticeably high for an object belong¬ 
ing to the B taxonomic class, which is generally supposed 
to include low-albedo a s teroid s. On the other hand, the 
I Shevchenko and Tedescol ll2006l) albedo of Pallas is also sub- 
stanti ally confirmed by m ore recent results based on WISE 
data (iMasiero et aklBoill ). The albedo value derived from 
the polarimetric slope turns out to be 0.088 ± 0.007, which 
would appear to be more in agreement with expectations for 
a body belonging to the low-albedo, B class. It is interesting 

^ If we call y = log(pv) and x = \og{h), by solving by a linear 
least-squares technique Eq.[2]and determining the corresponding 
error dy of y, it is easy to see that the corresponding error dpy 
turns out to be given by dpy = In(lO.O) py dy 


to note that [Pe Leon et al.l ll2012h observed a large sample 
of B-class objects and found a continuous variation of near- 
IR spectral slopes, possibly suggesting a variety of different 
compositions. This might be a consequence of the fact that 
the modern B class in cludes also some asteroids (the old 
F class of Tholen, see iTholen and Barucci I Il989l ) that in 
the 80s were kept separate based on their behaviour at the 
shortest wavelengths, which are no longer covered in the 
most modern spectroscopic investigations. One should also 
take into account that the surface properties of the largest 
asteroids like (2) Pallas, which retain a larger fraction of the 
material excavated in most impacts, can be different from 
those of smaller asteroids which lose a much larger frac¬ 
tion of the impact debris from most collisions to space. We 
already know, based on their different IR beaming param¬ 
eters, that the surfaces of large and smaller asteroids can 
be significantly different. In any case, polarimetry seems to 
indicate for (2) Pallas a low albedo. In the absence of any 
new, updated albedo and/or absolute magnitude value for 
this object coming from other sources, we accept this dis¬ 
crepancy. 

Another important problem which is most apparent in 
FiglSjconcerns the asteroid (64) Angelina, namely the object 
having t he highest albedo value i n our sample. The value 
given by IShevchenko and Tedescol ll2006l) for this object is 
0.474±0.047. According to our new calibration, the resulting 
albedo turns o ut to be 0.540 ± 0.085, margina lly in agree¬ 
ment with the [Shevchenko and Tedescj (l2006l ) value. The 
problem here is not this discrepancy per se, but rather the 
fact that, as shown in Fig.[2j in the high-albedo domain our 
new calibration tends to assign albedo values increasingly 
close to 1.0 to asteroids having increasingly shallower po¬ 
larimetric slopes. 

In Fig.[S]we show a comparison between our new slope- 
albedo relation (shown in red in the plot) and the most 
recent previous calibration, namely that of ICellino et al.l 
1I2OI2I) (black line). Though not being visually very different, 
an effect of the adopted log-log scale, the new calibration 
tends, at high values of albedo, to stay well above the val¬ 
ues predicted by the previous calibration. We will show in 
a separate paper that in the case of (44) Nysa, an asteroid 
belonging to the high-albedo E class, our new calibration 
assigns to this asteroid an albedo of about 0.9, which seems 
exceedingly high to be credible for a rocky body. In the 
domain of low-albedo asteroids, the larger steepness of the 
new best-fit line produces only a marginally better fit with 
respect to the previous calibration. 

This problem of an exceedingly high slope of our linear 
calibration seems to be connected with another major prob¬ 
lem that is apparent by looking at Fig. [21 This is the fact 
that at low albedo, the data look rather noisy, with some 
data points, having polarimetric slopes between 0.2%/deg 
and 0.3 %/deg, which are located well above or below the 
linear best-fit of the whole dataset. Trying to fit all the data, 
one is led to accept a linear best-fit whose steepness is a con¬ 
sequence of the presence of the lowest-albedo asteroids. This 
is a problem that has been known since a long time, and we 
will discuss it more extensively in the following Subsection. 
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Figure 6. The same as Fig.[2l but here the objects having albedo 
smaller than 0.08 (displayed here using open squares) were not 
used to derive the displayed linear best-fit. The corresponding 
values of the Ci and C 2 calibration coefficients are also indicated. 

4.1 Effects of excluding low-albedo objects 

In general terms, looking at Fig. [S] one could be tempted 
to conclude that it is hard to fit the whole dataset by using 
one unique linear relation. In particular, it may appear that, 
if one could drop a handful of objects having albedo lower 
than about 0.08, one could obtain a much better fit for the 
remaining objects. This is an old-debated subject, namely 
whether there is evidence of a saturation of the slope -albedo 
law at small albedo values. This kind of possible saturation 
is also much more evident in Pmin - albedo data, as we will 
see in Section [S] 

What happens if we exlude from our analysis low-albedo 
asteroids? Figure [6] shows the results of this exercise. As 
expected, the RMS deviation of the linear best-fit of the 
data, which is now much shallower than in the case in which 
we kept all the available measurements, is quite better, as 
shown in Table 4. One could conclude that excluding the 
asteroids having albedo lower than 0.08 is the best way to 
proceed to obtain an improved, and quite better, calibra¬ 
tion of the slope - albedo law. Figure [3 shows that in this 
way, the relative error on the derived albedo values turns 
out to be generally better than 20%, a quite good result. 
The problem, of course, is that in the practical applications 
of asteroid polarimetry, one has at disposal the polarimetric 
slope derived from observations, and wants to derive from 
it the albedo, which is unknown. Unfortunately, there is a 
range of values of polarization slope which is shared by ob¬ 
jects having albedo either around 0.05 or around 0.10. By 
using a slope - albedo relation which is not valid for low- 
albedo asteroids may produce a systematic overestimate of 
the albedo for dark objects. 

The best procedure to be adopted in practice may 


Figure 7. The same as Fig.|4l but here the objects having albedo 
smaller than 0.08 are not included in the analysis, and the plot 
shows the relative error of the albedo determinations. 
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Polarimetric Siope (%/deg) 


Figure 8. Results of a parabolic best-fit of all available slope - 
albedo data. The corresponding linear fits shown in Fig. [2] and in 
Fig.|6] are also displayed for a comparison (red and blue dashed 
lines, respectively). 
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be the following: when the polarimetric slope of an object 
is measured with good accuracy, it will be better to use 
the calibration coefficients obtained by dropping low-albedo 
objects, displayed in Table 1 and Fig. [6l but only when 
the polarimetric slope turns out to be smaller than about 
0.25%/deg. In this way, the relative error in the determina¬ 
tion of the albedo should be within 20%, a nice result, as 
shown in Fig. [T] 

If the slope is larger than the above value, and/or when 
the value of polarization slope is more uncertain, the best 
choice would be probably to use the calibration coefficients 
fitting the whole population, displayed in Table 1 and Fig. [2l 
The errors that one should expect for the higher values of 
polarimetric slopes should be in any case limited, of the 
order of about ±0.03, not negligible in relative terms, but 
in any case sufficient to correctly classify the objects as low- 
albedo asteroids. 

The polarimetric slope data shown in Fig. [2] can also 
suggest that a linear fit is not fully adequate to represent the 
whole dataset, and a parabolic fit could be more suited to 
better represent the data. This is shown in Fig. |8l in which 
a parabolic relation 

log(/i) = I{i(log(pv)f + H 2 log(pv) ± Hs 

is adopted, and the result of a best-fit procedure is shown. 
The linear plots already shown in Fig. [2] and Fig. |6] are also 
shown for a visual comparison. It is clear from the Figure 
that a parabolic fit applied to the whole dataset gives good 
results. The slightly worse fit of some objects having albedo 
around 0.2 (located below the parabolic fit) seems to be 
compensated by a much better fit of the polarimetric slope 
and albedo of the asteroid (2) Pallas, which is the most dis¬ 
crepant object when trying linear fits to represent the data. 
However, there is still the problem of the ambiguity affect¬ 
ing the values of albedo to be assigned to objects having 
polarimetric slopes between 0.2%/deg and 0.3%/deg, which 
cannot be solved. 

For this reason, we list the obtained values of the A, 
B, C coefficients of our parabolic fit in Table 1, but we 
are not claiming that there is such a big improvement with 
respect to the classical linear fit, specially when dropping 
the lowest-albedo objects, to force us to necessarily use a 
parabolic fit in the future. Things could change in the case 
that the discovery of other objects sharing the location of 
(2) Pallas in the h - pv plane would confirm that a linear 
fit is not really suited to adequately represent them, even 
by excluding from the analysis low-albedo objects. Another 
possibility is that future theoretical advances in the inter¬ 
pretation of light scattering phenomena could suggest that 
a parabolic fit is intrinsically more correct to fit the relation 
between polarimetric slope and geometric albedo, based on 
some physical arguments. 


5 MORE ON CLASSICAL APPROACHES: 
FITTING PHASE-POLARIZATION CURVES 

One can wonder whether the slope - albedo relation is really 
the best available choice to derive good estimates of aster¬ 
oid albedos. In fact, asteroid phase - polarization curves do 
not include only the (mostly) linear variation of Pi around 
the inversion angle. A negative polarization branch also ex- 


35 -1-1-1-1-■-■-■-1-^^^-1-■-1-1-[-■-■-■-1-^^-r 


30 


-a 


20 22 24 26 28 30 32 

A (%) 


Figure 9. The values of the best-fit parameters A and B found 
for the best 200 solutions of the linear-exponential relation ap¬ 
plied to polarimetric data of asteroid (216) Kleopatra, using a 
“genetic” algorithm. Note that all the plotted solutions for A and 
B give nearly identical residuals and produce best-fit curves char¬ 
acterized by values of h, Pmin other polarimetric parameters 
which are essentially identical (within the nominal error bars). 

ists, not to mention the behavior exhibited at large phase 
angles (not achievable for main belt objects) by near-Earth 
asteroids. 

Historically, another relation between albedo and polar¬ 
ization properties was found to involve the Pmin parameter, 
as shown in Eq. As mentioned in the previous Section, 
this has been generally abandoned in recent years, because 
Pmin data have been found by some authors to be quite scat¬ 
tered around the best-fit representation given by Eq.|3 On 
the other hand, one could also wonder whether this might 
be at least partly a consequence of the difficulty of deriving 
accurate values of Pmin from the observations, making use 
of visual extrapolations of rather sparse polarimetric data. 
This leads us to face the problem of finding suitable analyti¬ 
cal representations of the morphology of phase - polarization 
curves. 

In this paper we follow the example of previous authors, 
and we use the following exponential-linear relation: 

Pi = A(e~"l^ -l) + C-a (4) 

where a is the phase angle expressed in degrees, and A, 
B, C are parameters to be determined by means of best-fit 
techniques. This analytical representation has been found 
in the past to be suited to fit both phase - magnitude 
relations in asteroid photometry, and phase - polariza¬ 
tion curves in asteroid polarimetry (iKaasalainen et ^l2003l : 
iMuinonen et al.ll2009lf . Some examples of practical appli¬ 
cations of the above relation are the best-fit curves of the 
phase-polarization curves of the asteroids shown in Fig. [T] 
According to its mathematical representation, when the pa- 


© 2014 RAS, MNRAS OOO.ITIfTSl 







10 


A. Cellino et al. 


■o 




o 


0.75 


0.7 


0.65 


0.6 


0.55 


0.5 



20 22 24 26 28 30 32 

A (%) 


Figure 10. The same as in Fig.O but for the A and C parameters 
in Eq. |4] 


rameters A, B, C are all positive, the exponential-linear re¬ 
lation describes a curve characterized by a negative polar¬ 
ization branch between 0 and an inversion angle Qinv The 
trend tends to become essentially linear at large phase an¬ 
gles, where the exponential term tends quickly to zero. 

The computation of the best-fit representation of any 
phase-polarization curve using Eq. |T] can be done in many 
ways. In the present paper, we use a genetic algorithm, 
which, starting from a random set of A, B, C values, explores 
the space of possible solution parameters and finds the set 
of A, B, C values producing the smallest possible residuals. 
Due to the intrinsic properties of a genetic approach, the 
algorithm is launched several times, in order to have a cor¬ 
respondingly high number of solutions, in order to ensure 
that we are not missing the best possible one. 

It should be noted, however, that the evaluation of the 
errors of any polarimetric parameter derived by a best-fit of 
the phase - polarization curve using Eq.jT] is complicated by 
the fact that the parameters A, B and C that minimise the 
are correlated. To illustrate this situation. Fig. [5] shows 
that different pairs of A and B values produce fits nearly 
indistinguishable in terms of r.m.s. residuals (the differences 
being not larger than 0.0015). Figure [in] shows a similar sit¬ 
uation for the parameters A and C. This means t hat the non 
diago nal elements of the error matrix (see, e.g.. ISevingtonl 
Il969l i are not negligible with respect to the diagonal ones, 
and therefore the calculation of the error on any polarimetric 
parameter derived from the exponential-linear curve should 
take into account all the various co-variances. Our method 
for the minimization is based on a genetic algorithm, 
which does not produce automatically the error matrix. For 
the estimate of the error on polarimetric parameters derived 
by best-fit values of A, B, C, it is therefore more practical to 
adopt an alternative method. Let us assume, as an example. 


that we are interested in determining Pmin and its corre¬ 
sponding uncertainty. The resulting value of Pmin will be 
the one obtained using the A, B, C values giving the small¬ 
est x^. As for the error to be assigned to this determination 
of Pmin, the method that we adopt consists of calculating all 
Pmin values corresponding to identified sets of A,B,C pa¬ 
rameters that produce a fit of the phase - polarization curve 
giving x^ values such that ^ Xmin + 1- We then define 
as error on Pmin the half difference between the extremes 
of the various Pmin values so obtain ed. This approach was 
followed e.g. bv iBagnulo et all || l995|) and i s cons istent with 
the error analysis presented bv iBevingtonI (Il969li . Needless 
to say, this procedure may be applied to the determination 
of any polarimetric parameter (other than Pmin) derived by 
an exponential-linear fit of the phase-polarization curve. 

In what follows, since we want to limit our analysis 
only to high-quality and well-covered phase - polarization 
curves, we impose some strict constraints on the selection 
of the objects for which we compute a best-fit using the 
exponential-linear relation. In particular: 

• we exclude a priori from our analysis all measurements 
having a nominal accuracy of Pi worse than 0.20. 

• we require to have at least 4 accepted measurements 
taken at phase angles > 2°. 

• we require to have at least 1 accepted measurement 
taken at phase angles > 17°. 

• we require to have at least 1 accepted measurement 
taken at phase angles < 14°. 

• we require to have at least 3 accepted measurements 
taken at phase angles < 30°. 

As in the case of the computation of the polarimetric 
slope described in the previous Section, we limit our analy¬ 
sis to available data obtained in V hlter. Here, however, we 
add an additional constraint: we do not use for calibration 
purposes the best-fit phase - polarization curves of asteroids 
for which we have fewer than 10 accepted Pi measurements. 
All the criteria described above are dictated by our will to 
restrict our calibration procedures to the objects having ex¬ 
cellently determined and optimally sampled phase - polar¬ 
ization curves, only. 


5.1 Use of Pmin 

Having at disposal a best-fit representation of a phase- 
polarization curve according to Eq. [I] one can compute 
the resulting Pmin value and the corresponding phase an¬ 
gle a(Pmin) at which it is found. More in detail, a(Pmin) can 
be computed by equaling to zero the first derivative of Eq.|4l 
from which we obtain: 

a(Pmin) = 

then, Pmin can be computed as Pr(a(Pmin)) using Eq.|Tl The 
nominal uncertainty in Pmin has to be computed by doing a 
formal propagation of the errors of the A, B, C parameters 
found in the best-fit solution of xEq. [T] using the procedure 
explained above. 

In the case of the Pmin - albedo relation described by 
Eq. [3] the results of our exercise are shown in Fig . 1111 in¬ 
cluding the best-fit values we find for the Cs and C 4 coeffi¬ 
cients. The calibration is based on data of the Pmin values of 
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Figure 11. Best-fit relation between Pmin and geometric albedo 
for an available sample of 20 asteroid targets included in the 
IShevchenko and Tedesc3 ll2006l) list. Only 16 objects for which 
we have at disposal at least ten polarimetric measurements have 
been used to derive the best-fit solution. These objects are plot¬ 
ted as full, black symbols. Open, green symbols represent objects 
for which we have fewer than ten observations, and were not used 
in the least-squares computation. 


16 objects belonging to the lShevchenko and Tedescol (l2006l l 
list, for which we have at least 10 polarimetric measure¬ 
ments. Four additional asteroids, for which we have fewer 
than ten measurements, are also shown using different sym¬ 
bols, but they were not used in the computation of the 
best-fit. The obtained values of C 3 and Ci calibration co¬ 
efficients, together with their corresponding errors, are also 
given in Table 1. We see that, not unexpectedly, the dis¬ 
tribution of the points in the Pmin - albedo plane, in our 
log-log plot (in Fig lllll makes it difficult to find a satisfac¬ 
tory linear fit. Correspondingly, the agreement of the re¬ 
sultin g albedos with those of the IShevchenko and Tedesca 
(|200(1) list is significantly worse than in the case of the cali¬ 
bration based on the polarimetric slope. As shown in Table 
[6l we find a large discrepancy in the case of the bright as¬ 
teroid (64) Angelina, for which a very high value of albedo 
of 0.600 ± 0.044 is obtained from the Pmin - albedo rela- 
tion, much larger than t he 0. 474 ± 0.047 value listed by 
IShevchenko and Tedescol (1200(1 ). The majority of the other 
asteroids of our sample, conversely, tends to have an albedo 
underestimated with respect to the calibration values, apart 
from a few ones having the lowest albedo values. Among 
them, we note the difference between the very low albedo 
value found by IShevchenko and Tedescol (1200(1) for asteroid 
(444) Gyptis (0.037 ± 0.004) and the value of 0.106 ± 0.002 
that we obtain from its Pmin value. This object, however, was 
not used in the derivation of the best-fit computation, be¬ 


Figure 12. The same as Fig. 1111 but here the objects having 
albedo lower than 0.08 (sowed using open symbols) were not used 
in the computation of the linear best-fit. 

cause only 6 polarimetric measurements are currently avail¬ 
able for it, and moreover the phase - polarization curve (not 
shown) is quite noisy. We note also that the albedo value 
found by IShevchenko and Tedescol (1200(1) is very low, and 
might possibly be slightly underestimated. 

By simply looking at Fig llll a saturation of the Pmin 
- albedo relation at low albedo values is even more evident 
than in the case of the slope - albedo relation analyzed in 
Section KB It seems therefore that a removal of asteroids 
having albedo less than about 0.08 from the best-ht compu¬ 
tation is even more justihed in this case. The result of this 
exercise is shown in Fig. 1121 The improvement of the resid¬ 
uals, listed in Table 4, is very important, as also visually 
shown in the Figure. 

The big improvement of the obtained best-ht makes 
this Pmin - albedo relation much more suitable for the de¬ 
termination of the albedo, but, again, this refers to only a 
more limited interval of possible Pmin values, in particular 
those lower (in absolute value) than about 1%. For objects 
having deeper Pmin, the corresponding interval of possible 
albedo values is exceedingly wide to be used to derive a use¬ 
ful albedo determination. Based on our results, we conhrm 
therefore that, in general, the use of Pmin as a reliable diag¬ 
nostic of the albedo, but for asteroids exhibiting a shallow 
polarization branch, should not be encouraged. 

5.2 An alternative derivation of the polarimetric 
slope 

Using the global htting of the phase polarization curves 
given by Eq.|4l it is also possible to modify the way to derive 
the polarimetric slope h. In so doing, one can make use of 
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Figure 13. The same as Fig. [J] but here the polarimetric 
slope hABC is computed as the first derivative of Pr (accord- 
ing to Eg.ro at the inve r sion a ngle otin^ for 20 asteroids of the 
IShevchenko and Tedesc'3 ll2nnol 'l list. Four objects having fewer 
than ten polarimetric measurements were not used for the com¬ 
putation of the best-fit. They are indicated by open, green sym¬ 
bols. 

all the available polarimetric measurements, and not only of 
those obtained in a more or less narrow interval of phase an¬ 
gles centered around the inversion angle. In particular, the 
polarimetric slope can be determined as the first derivative 
of Pr with respect to the phase angle (using Eqj4)| , where the 
derivative has to be computed at the inversion angle ainv. 
In turn, the value of ainv can also be derived with excellent 
accuracy from Eq.|4l We adopt here a very simple numeri¬ 
cal approach. Having determined the values of parameters 
A, B, C corresponding to the nominal best-fit solution, we 
make an iterative computation of Pr starting from an initial 
phase angle value of 1 degree, and using a fixed increment of 
-1-0.02 deg in phase. When Pr(i-|- 1) •Pr(*) becomes negative, 
we consider that ainv is equal to Pr(i) + 0.01 degrees, with 
an uncertainty of 0.02 degrees. 

Having determined the inversion angle, we can com¬ 
pute the resulting polarimetric slope as the first derivative 
of Eq.|4l using the same procedure already adopted for Pmin 
to derive its nominal uncertainty. In what follows we will 
always call hABC these new value of the slope computed 
as explain ed above. We obtained hABC v alues for 20 ob¬ 
jects of the lShevchenko and Tedescol ll2006li list, as shown in 
Fig |13l Also in this case, however, we did not use in the best- 
fit computation the data of four objects having fewer than 
10 polarimetric measurements. These four asteroids, namely 
(78) Diana, (216) Kleopatra, (431) Nephele and (444) Gyp- 
tis, are indicated in Fig |13l bv means of open, green symbols. 
Note that our sample includes now three extra objects, (27) 
Euterpe, (41) Daphne and (47) Aglaja, which did not sat¬ 
isfy our previous acceptability criterion for the computation 


Figure 14. The same as Fig. 1131 but here objects having an 
albedo lower than 0.08 (shown using open squares) are not used 
in the computation of the best-fit. 


of the polarimetric slope h carried out in Section [4] Con¬ 
versely, in Section |4] we made use of polarimetric slopes for 
the two asteroids (105) Artemis and (124) Alkeste, which do 
not satisfy our criteria for the computation of Habc- 

The resulting values of the calibration coefficients Ci 
and C 2 are shown in Fig. 1131 and they are also listed, to¬ 
gether with their errors, in Table 1. As in the cases seen 
above, we also computed the best-fit values of the calibra¬ 
tion parameters which are obtained by removing from the 
computation the asteroids having albedo lower than 0.08. 
The results of this exercise, listed in Table 1, are also shown 
in Fig. [141 

By looking at the results, we find that the linear best-fit 
of the slope - albedo relation look, again, reasonably good. 
However, there is not any improvement with respect to the 
case when the polarimetric slope was computed by doing 
a more trivial linear fit of the available data around the 
inversion angle (see Fig.[3l). As opposite, the RMS d eviation 
with respect to the lShevchenko and Tedescol (l2006li albedos 
turn out to be slightly worse, as shown in Table 4. 

The computation of the polarimetric slope from a sim¬ 
ple linear fit of data distributed around the inversion an¬ 
gle, or from the computation of the first derivative of 
PrComputed at the inversion angle, which would better cor¬ 
respond to the ideal definition of this parameter, is therefore 
not fully equivalent. It turns out that, opposite to our own 
expectations, the simpler (purely linear) approach seems to 
give slightly better results, in spite of all the uncertainties. 

The results of the exercises described in this Section are 
summarized in Table 2, in which we list all the polarimet¬ 
ric parameters considered in our analysis (including some 
which will be exp lained in the next Section) , and i n Table 3, 
where we list the IShevchenko and Tedescol ll2006h values of 
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Table 1. Resulting values and nominal uncertainties of the calibration coefficients in different albedo - polarization relations considered 
in this paper. The result of a possible parabolic fit of the polarimetric slope as a function of albedo is also listed (see text). 


log(pv) = Cl log(h) + C2 

log^v-) = Cl log(h) + C2 (pv Ss 0 . 08 ) 

log(Pv) = Cl log(/iAsc) + C2 

log(pv) = Cl logihABc) + C2 (pv 5 = 0 . 08 ) 

log(pv) = C3 log(Pniin) + C4 

log(pv) = C'a log(Pmin) + C4, (pv ^ 0 . 08 ) 

log(pv) = C^i log( 5 ') + C^2 

log(pi/) = C *1 p* + C*2 

log(/i) = Hi(log{pv)f + H2 log(pv) + P3 


Cl = -1.111 ± 0.031 
Cl = -0.800 ± 0.041 
Cl = -1.139 ± 0.026 
Cl = -0.780 ± 0.037 
C 3 = -1.419 ± 0.034 
C 3 = -0.869 ± 0.042 
C^i = -0.987 ± 0.022 
C*i = -0.896 ± 0.029 
Hi = -1.294 ±0.001 


C 2 = -1.781 ±0.025 
C 2 = -1.467 ±0.037 
C 2 = -1.850 ±0.021 
C 2 = -1.469 ±0.036 
C 4 = -0.918 ± 0.006 
Ca = -0.789 ± 0.008 
C^2 = -0.458 ± 0.013 
C*2 = -1.457 ±0.018 

P 2 =-3.140 ±0.001 P 3 =-2.428 ± 0.001 



Fi gure 15. The p*-albedo relation for the 13 asteroids of the list 
of Ishevchenko and Tedescol ll2006l ) for which a reliable estimate 
of p* is possible based on obtained values of the h polarimetric 
slope and of Pmin- The best-fit solution is plotted together with 
the individual data. The newly derived values of the A and B 
parameters and their resulting errors are also indicated. 


albedo, together with the corresponding values of albedo de¬ 
rived from the considered polarimetric parameters and their 
nominal errors. Note that in Table 2, we always give for the 
inversion angle ainv the value obtained from the best-fit of 
the whole phase - polarization curve using the exponential - 
linear relation. Only in cases in which this is not available, 
but we have at disposal a polarimetric slope h obtained as 
described in Section [4] we assign to ainv the value corre¬ 
sponding to the intersection of the polarimetric slope with 
the Pt — 0 line. 

The polarimetric parameters h, Habc and Pmin, and 
corresponding albedos obtained from calibrations based on 
all and only the objects having albedos larger than 0.08, as 
discussed in this and the previous Section, are also listed 
in Table 5. The improvement of the agreement between the 
albedos obtained from polarimetric parameters and the albe- 


Figure 16. The H'-a lbedo relation in log - lo g scale for 20 as¬ 
teroids of the list of Ishevchenko and Tedes^ (l2006l l for which 
we were able to compute the 11' parameter. Four asteroids hav¬ 
ing fewer than 10 polarimetric measurements were not used in 
the computation of the linear best-fit and are indicated by open, 
green symbols. The best-fit solution corresponding to the relation 
log(pv) = Cl log(\I') ± C 2 is plotted together with the individ¬ 
ual data. The best-fit values of the parameters and their nominal 
uncertainties are also indicated in Table 1. 

dos given bv I Shevchenko and Tedescol (I 2 OO 6 I) is evident. Our 
considerations about the best possible use of these calibra¬ 
tions have been already exposed in previous Sections. 


6 OTHER POLARIMETRIC PARAMETERS 

The failure of our attempt to obtain more accurate albedo 
values by using a polarimetric slope (/iasc) obtained by a 
formal computation of the first derivative of Eq. |4] at the 
inversion angle, can be important. A much simpler linear 
fit of polarimetric data spread over a large interval of phase 
angles, seems to be capable of giving slightly more accurate 
albedo solutions. This can be an indication that using po- 
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Py (Shevchenko & Tedesco 2006) 


Figure 17. Differences between the albedos of 20 objects, derived 
using our calibration of the \E' - albedo relation, and the corre - 
sponding albedo values given bv IShevchenko and Tedescol ll2006tl . 
Four asteroids having fewer than 10 polarimetric measurements 
were not used in the computation of the linear best-fit and are 
indicated by open, green symbols. 


larimetric parameters describing the behavior of the phase - 
polarization curve only at some single value of phase angle, 
or in a limited portion of the phase angle interval, for the 
determination of the geometric albedo, could be not a very 
good idea. 

According to current evidence, we know that low albedo 
asteroids exhibit a deeper value of Pmin as well as a steeper 
linear polarization slope over a large interval of phase an¬ 
gles. Asteroids having increasingly higher albedos exhibit 
an opposite behavior (increasingly shallower Pmin and gen¬ 
tler h). There are also differences in the typical values of 
the inversion angle for different classes of objects, as found, 
as an example , in t he case of the F taxonomic class by 
iBelskava et al.l ll2005ll . One can imagine many different ways 
to attempt a new calibration of the albedo - polarization re¬ 
lationship trying to exploit the above evidence and make 
use of the overall morphology of the phase - polarization 
curve. A first attempt in this direction was proposed by 
iMasiero et al.l (120121 '). 

These authors proposed to use a new observable, they 
called p* (p-star), defined as a parameter of maximum po¬ 
larimetric variation, given by 

p* = Wl log(/l) -I- W2 log(Pmin) 

where h is, again, the classical polarimetric slope, and Wi 
and W 2 are two pa rameters whose values were found by 
IMasiero et al.l (|2ni2l l to be Wi = 0.79 ± 0.02 and W 2 = 
0.61 ± 0.03. The authors used in their analysis a data-set of 
177 asteroids having an albedo value estimated from thermal 
radiometry data produced by the WISE mission, whereas 
polarimetric data were taken by the authors from the liter¬ 


ature. For 65 asteroids of this sample, the authors derived 
a value of h, whil e for 112 objec t s the y derived a value of 
Pmin. In this way, IMasiero et al.l (l2012h derived a new cal¬ 
ibration of the albedo - polarization relationship based on 
WISE albedos and the newly introduced p* parameter. The 
relation they found was: 

log(pv) = C *1 p* +C*2 


and they found for the C*i and C*2 coefficients the values 
C*i = -1.04 ± 0.04 and C*2 = -1.58 ± 0.09. 

Since we use in our analysis a much smaller number 
of asteroids having presumably more accurate albedo values 
not derived from thermal radiometry, and we use a different 
set of polarimetric data including a large number of previ¬ 
ously unpublished observations, we have decided to derive a 
new calibration of the albedo - p* relation, using the data 
at our disposal, while keeping the same defi nition of the p* 
parameter as given by IMasiero et al.l (l2012f). In partic ular, 
we kept the values computed by Masiero et al.l (l2012f) for 
the Wl and W2 parameters, which were derived by the au¬ 
thors based on their own analysis of polarimetric data. We 
computed then the p* parameter for our sample of asteroids 
using Pmin values already obtained from our best-fit of Eq.|4] 
and polarimetric slopes h derived as described in Section [T] 
We had at disposal estimates for both h and Pmin for 13 ob¬ 
jects, only. The results of this exercise are shown in Fig. 1151 
in which we also show the new values that we found for the 
C*i and C*2 parameters. In particular, we have: 


C*i = -0.896 ±0.029 


and 


C*2 = -1.457 ±0.018 

The resulting fit is fairly good, and we confirm that p* is an¬ 
other useful parameterfor albedo determination. The result¬ 
ing RMS deviations, however, are higher than those corre¬ 
sponding to the slope - albedo relation discussed above, both 
using either h or habc, as shown in Table 4. We find a quite 
big discrepancy concerning the p redicted albedo for (64) 
Angel ina and the corresponding IShevchenko and Tedesca 
J2006tl value. In the region of higher p* values (right region 
of the plot), moreover, the scatter of the albedos around the 
best-fit line is fairly high. 

There are, of course, other possibilities to use polari¬ 
metric parameters describing the overall morphology of the 
observed phase - polarization curves. A very simple idea is 
to use some parameter built directly from the obtained val¬ 
ues of polarization Pr taken at very different values of phase 
angle, possibly including both the negative and the positive 
polarization branches. In this paper, we introduce such a 
parameter, that we call tk, and we define it as 

4' = Pr(30°)-Pr(10°) 

where the dependence of Pr upon the phase angle is assumed 
to be given by Eq.[4l We derived t herefore the T parameter 
for ou r sample of asteroids in the I Shevchenko and Tedescol 
(I2OO6II list following the same criteria already described 
above in the case of our computations of Pmin and hABC- 
This corresponds to 20 asteroids, listed in Table 2. For the 
purposes of calibration, again, we did not use data of four 
asteroids having fewer than 10 polarimetric measurements, 
which in Fig. [1^ are displayed using different symbols. The 
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Table 2. Summary of the polarimetric parameters found for all asteroids included in the lShevchenko and Tedescol l|2006l) list, for which 
we have a suitable coverage of the phase - polarization curves. Each asteroid is identified by its number N. The second column gives the 


number A^obs of polarimetric measurements used in the analysis. For the meaning of the other parameters, see the text. 


N 

-^obs 


h 


a 

inv 


o(^’min) 

P 

min 


4' 


hABC 


pu 


1 

33 

0.2549 

± 

0.0041 

18.13 

± 

0.02 

7.20 

± 

0.02 

-1.683 

± 

0.005 

4.863 

± 

0.020 

0.2467 

± 

0.0009 

0.331 

± 

0.015 

2 

22 

0.2223 

± 

0.0015 

18.57 

± 

0.02 

7.59 

± 

0.07 

-1.413 

± 

0.017 

4.123 

± 

0.010 

0.2111 

± 

0.0004 

-0.424 

± 

0.014 

3 

26 

0.1029 

± 

0.0039 

20.31 

± 

0.02 

8.04 

± 

0.08 

-0.732 

± 

0.009 

1.724 

± 

0.023 

0.0949 

± 

0.0013 

-0.863 

± 

0.024 

4 

26 

0.0691 

± 

0.0012 

22.37 

± 

0.02 

9.22 

± 

0.08 

-0.558 

± 

0.005 

1.150 

± 

0.006 

0.0701 

± 

0.0005 

-1.071 

± 

0.025 

8 

28 

0.1067 

± 

0.0016 

20.05 

± 

0.02 

8.30 

± 

0.03 

-0.671 

± 

0.007 

1.743 

± 

0.006 

0.0950 

± 

0.0003 

-0.873 

± 

0.021 

27 

12 


- 


21.51 

± 

0.02 

7.27 

± 

0.22 

-0.593 

± 

0.010 

1.062 

± 

0.033 

0.0582 

± 

0.0023 


- 


39 

20 

0.0745 

± 

0.0093 

21.43 

± 

0.02 

8.66 

± 

0.07 

-0.712 

± 

0.013 

1.563 

± 

0.017 

0.0906 

± 

0.0009 

-0.981 

± 

0.049 

41 

11 


- 


22.53 

± 

0.02 

10.54 

± 

0.04 

-1.664 

± 

0.015 

4.040 

± 

0.034 

0.2591 

± 

0.0019 


- 


47 

11 


- 


17.95 

± 

0.02 

7.90 

± 

0.20 

-1.378 

± 

0.031 

5.037 

± 

0.312 

0.2412 

± 

0.0109 


- 


51 

20 

0.2713 

± 

0.0040 

20.35 

± 

0.02 

8.29 

± 

0.05 

-1.950 

± 

0.015 

4.784 

± 

0.023 

0.2644 

± 

0.0013 

-0.271 

± 

0.015 

64 

13 

0.0434 

± 

0.0043 

18.77 

± 

0.02 

6.75 

± 

0.26 

-0.323 

± 

0.014 

0.757 

± 

0.028 

0.0395 

± 

0.0014 

-1.376 

± 

0.048 

78 

5 


- 


22.19 

± 

0.02 

10.42 

± 

0.29 

-1.497 

± 

0.048 

3.835 

± 

0.481 

0.2394 

± 

0.0307 


- 


85 

11 

0.2648 

± 

0.0071 

19.07 

± 

0.02 

8.68 

± 

0.16 

-1.375 

± 

0.017 

4.825 

± 

0.108 

0.2414 

± 

0.0034 

-0.372 

± 

0.016 

105 

6 

0.2872 

± 

0.0032 

19.75 

± 

0.31 


- 



- 



- 



- 



- 


124 

5 

0.0745 

± 

0.0041 

19.79 

± 

1.46 


- 



- 



- 



- 



- 


129 

10 

0.0953 

± 

0.0077 

21.07 

± 

0.02 

7.61 

± 

0.03 

-0.849 

± 

0.011 

1.676 

± 

0.009 

0.0928 

± 

0.0004 

-0.850 

± 

0.035 

216 

9 


- 


18.83 

± 

0.02 

8.93 

± 

0.09 

-1.028 

± 

0.011 

4.151 

± 

0.124 

0.1971 

± 

0.0041 


- 


230 

15 

0.1045 

± 

0.0047 

20.45 

± 

0.02 

7.55 

± 

0.06 

-0.937 

± 

0.012 

1.991 

± 

0.017 

0.1089 

± 

0.0010 

-0.792 

± 

0.025 

324 

20 

0.2890 

± 

0.0046 

19.73 

± 

0.02 

8.70 

± 

0.05 

-1.655 

± 

0.016 

4.987 

± 

0.024 

0.2644 

± 

0.0010 

-0.292 

± 

0.014 

431 

7 


- 


19.67 

± 

0.02 

9.29 

± 

0.07 

-1.360 

± 

0.024 

4.864 

± 

0.048 

0.2478 

± 

0.0023 


- 


444 

6 


- 


20.65 

± 

0.02 

9.77 

± 

0.04 

-1.098 

± 

0.012 

3.477 

± 

0.041 

0.1914 

± 

0.0021 


- 


704 

32 

0.3074 

± 

0.0054 

15.73 

± 

0.02 

7.02 

± 

0.01 

-1.310 

± 

0.006 

6.478 

± 

0.011 

0.2692 

± 

0.0006 

-0.333 

± 

0.012 


Table 3. Resulting albedo values py for all asteroids belonging to the lShevchenko and Tedescol 1I2OO6II list, for which we have polarimetric 
observations suitable to derive an albedo value using one or more of the techniques explained in the text. Only the albedo corresponding 
to the calibration co mputed using all available objec ts (not only those with py > 0.08) are listed. The last column gives, for a comparison, 
the albedo given bv IShevchenko and Tedescol l|2006h . 


Number Pv{h) Pv{hABc) PviPmin) Pv{'^) Pv{p^) py(S(^T) 


1 

0.076 

± 

0.006 

0.070 

± 

0.004 

0.058 

± 

0.001 

0.073 

± 

0.003 

0.069 

± 

0.004 

0.094 

± 

0.006 

2 

0.088 

± 

0.007 

0.083 

± 

0.005 

0.074 

± 

0.002 

0.086 

± 

0.004 

0.084 

± 

0.005 

0.145 

± 

0.009 

3 

0.207 

± 

0.021 

0.206 

± 

0.016 

0.188 

± 

0.005 

0.203 

± 

0.007 

0.207 

± 

0.018 

0.187 

± 

0.006 

4 

0.322 

± 

0.033 

0.292 

± 

0.025 

0.276 

± 

0.008 

0.303 

± 

0.009 

0.318 

± 

0.031 

0.350 

± 

0.020 

8 

0.199 

± 

0.018 

0.206 

± 

0.016 

0.213 

± 

0.005 

0.201 

± 

0.007 

0.212 

± 

0.018 

0.197 

± 

0.012 

27 


- 


0.360 

± 

0.036 

0.254 

± 

0.008 

0.328 

± 

0.014 


- 


0.298 

± 

0.045 

39 

0.297 

± 

0.051 

0.218 

± 

0.017 

0.196 

± 

0.006 

0.224 

± 

0.007 

0.264 

± 

0.034 

0.246 

± 

0.015 

41 


- 


0.066 

± 

0.004 

0.059 

± 

0.001 

0.088 

± 

0.004 


- 


0.061 

± 

0.004 

47 


- 


0.071 

± 

0.006 

0.077 

± 

0.003 

0.071 

± 

0.005 


- 


0.060 

± 

0.002 

51 

0.071 

± 

0.005 

0.064 

± 

0.004 

0.047 

± 

0.001 

0.074 

± 

0.003 

0.061 

± 

0.003 

0.097 

± 

0.006 

64 

0.540 

± 

0.085 

0.560 

± 

0.059 

0.600 

± 

0.044 

0.458 

± 

0.022 

0.597 

± 

0.084 

0.474 

± 

0.047 

78 


- 


0.072 

± 

0.011 

0.068 

± 

0.003 

0.092 

± 

0.012 


- 


0.086 

± 

0.003 

85 

0.072 

± 

0.006 

0.071 

± 

0.004 

0.077 

± 

0.002 

0.074 

± 

0.004 

0.075 

± 

0.004 

0.054 

± 

0.003 

105 

0.066 

± 

0.005 


- 



- 



- 



- 


0.047 

± 

0.005 

124 

0.297 

± 

0.034 


- 



- 



- 



- 


0.240 

± 

0.036 

129 

0.226 

± 

0.029 

0.212 

± 

0.017 

0.152 

± 

0.004 

0.209 

± 

0.007 

0.202 

± 

0.020 

0.183 

± 

0.018 

216 


- 


0.090 

± 

0.006 

0.116 

± 

0.002 

0.085 

± 

0.004 


- 


0.170 

± 

0.010 

230 

0.204 

± 

0.021 

0.177 

± 

0.013 

0.132 

± 

0.003 

0.177 

± 

0.006 

0.179 

± 

0.015 

0.192 

± 

0.019 

324 

0.066 

± 

0.005 

0.064 

± 

0.004 

0.059 

± 

0.002 

0.071 

± 

0.003 

0.064 

± 

0.003 

0.051 

± 

0.003 

431 


- 


0.069 

± 

0.004 

0.078 

± 

0.002 

0.073 

± 

0.003 


- 


0.122 

± 

0.018 

444 


- 


0.093 

± 

0.006 

0.106 

± 

0.002 

0.102 

± 

0.004 


- 


0.037 

± 

0.004 

704 

0.061 

± 

0.004 

0.063 

± 

0.004 

0.082 

± 

0.001 

0.055 

± 

0.003 

0.069 

± 

0.004 

0.057 

± 

0.002 
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Table 4. Average RMS deviation of polarimetrically-derived albedo with respect to the values in the IShevchenko and Tedescol i|20Q6l ') 
list, using different possible albedo - polarization relations described in the text. For each case, N is the number of observed asteroids 
used to obtain the calibration 


Albedo computed from: 


RMS deviation 

N 

h slope from linear fit 

(all objects) 

0.035 

15 

h slope from linear fit 

(only objects having py > 0.08) 

0.033 

11 

habc slope from linear-exponential fit 

(all objects) 

0.038 

16 

habc slope from linear-exponential fit 

(only objects having py > 0.08) 

0.034 

11 

P ■ 

min 

(all objects) 

0.051 

16 

P . 
min 

(only objects having py > 0.08) 

0.035 

11 

$ = Pr(30) - Pr(lO) 

(all objects) 

0.026 

16 

P* 

(all objects) 

0.043 

13 


Table 5. List of polarimetric parameters h, Habc ^nd Pmin obtained for 15 asteroids of the lShevchenko and Tedescol l|2006l ') list, having 
albedo larger than 0.08, and the corresponding albedos computed using the calibrations based on these objects, only (see Table 1). Ntot 
is the number of polarimetric o bservations available for each object. The last column gives, for a comparison, the albedo value listed by 
IShevchenko and Tedescol ||2006|) . 


N Ntot h 


hABC ^min 


pvih) 


Pv(hABc) Pv(f’min) Pv(ST) 


1 

33 

0.2549 ± 0.0041 

0.2467 ± 0.0009 

-1.683 ±0.005 

2 

22 

0.2223 ±0.0015 

0.2111 ± 0.0004 

-1.413 ±0.017 

3 

26 

0.1029 ± 0.0039 

0.0949 ± 0.0013 

-0.732 ±0.009 

4 

26 

0.0691 ± 0.0012 

0.0701 ± 0.0005 

-0.558 ±0.005 

8 

28 

0.1067 ±0.0016 

0.0950 ± 0.0003 

-0.671 ±0.007 

27 

12 

- 

0.0582 ± 0.0023 

-0.593 ±0.010 

39 

20 

0.0745 ± 0.0093 

0.0906 ± 0.0009 

-0.712 ±0.013 

51 

20 

0.2713 ± 0.0040 

0.2644 ± 0.0013 

-1.950 ±0.015 

64 

13 

0.0434 ± 0.0043 

0.0395 ± 0.0014 

-0.323 ±0.014 

78 

5 

- 

0.2394 ± 0.0307 

-1.497 ±0.048 

124 

5 

0.0745 ± 0.0041 

- 

- 

129 

10 

0.0953 ± 0.0077 

0.0928 ± 0.0004 

-0.849 ±0.011 

216 

9 

- 

0.1971 ± 0.0041 

-1.028 ±0.011 

230 

15 

0.1045 ± 0.0047 

0.1089 ± 0.0010 

-0.937 ±0.012 

431 

7 

- 

0.2478 ± 0.0023 

-1.360 ±0.024 

Figure 

shows 

a log - log: plot 

of d' as a function of the 


IShevchenko and Tedescol (I 2 OO 6 I I albedo, according to the re¬ 
lation 

log(pv) = C^i log(4') -I- C ^2 

identical to the classical slope - albedo relation, but using 
the new parameter instead of the polarimetric slope h. 
Fig. 1161 shows the best-fit linear solution, and the resulting 
values of the Cipi and Cip 2 parameters, which are also listed, 
together with their corresponding uncertainties, in Table 1. 

The fit looks good, and this is confirmed by the result¬ 
ing RMS deviations, which are found to be slightly lower 
than in the case of all other albedo - polarization rela¬ 
tions examined in this paper, including those obtained by 
excluding low-albedo objects, as shown by the in Table 4. 
In Fig. 1171 we plot the differences between the albedo val¬ 
ues prod uced by our 4^ - albedo ca li bratio n and the albedos 
given by IShevchenko and Tedescol 1 I 2 OO 6 II . The albedo val¬ 
ues obtained from the 41 p arameter tend to be always very 
close to the corresponding IShevchenko and Tedescol (l2006l i 
values, although the relative uncertainty tends to be fairly 
high for the darkest asteroids. 

With respect to the results obtained by using the polari¬ 
metric slope as shown in previous Sections, it is encouraging 
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to note that our 41 - albedo rel ation seems to fit nicely the 
IShevchenko and Tedescol (l2006[l albedos in the whole inter¬ 
val covered by the data. The albedo obtained for (64) An¬ 
gelina seems to suggest that we have no longer the problem 
of a possible overestimation of the albedo of bright objects, 
which affected the calibration of the classical slope - albedo 
law, as seen in Section [d] We also note that the value of 
Cipi = —0.987 ± 0.022, obtained in our calibration of the 
41 - albedo relation, is also formally in agreement, within 
the uncertainty, with an even simpler hyperbolic relation 
pv = A'/4', with K = 10*^^ ~ 0.348. 

We find, again, some problems with asteroid (2) Pal¬ 
las, for which an albedo of 0.086 ± 0.004 is found by our 
4^ calibr ation, still much lower th a n the 0.145 albedo indi¬ 
cated by IShevchenko and Tedescol 1 I 2 OO 6 II . The same is true 
for (216) Kleopatra, which has a 4* value nearly identical to 
that of Pallas, corres ponding to an albedo of 0 . 085 ± 0.004, 
much lower than the lShevchenko and Tedescol ll2006li value 
of 0.170±0.010. As opposite, for (444) Gyptis we have prob¬ 
lems of a remarkable overestimation of th e albedo, but in this 
case the lShevchenko and Tedescol (l2006l i albedo is extremely 
low, 0.037 ±0.004. Both Kleopatra and Gyptis are asteroids 
for which we have fewer than 10 polarimetric measurements. 
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In the case of (444) Gyptis, the phase - polarization curve 
is quite noisy, whereas this is not the case for Kleopatra. 


7 CONCLUSIONS AND FUTURE WORK 

The solution of the problem of determining the best possible 
calibration of the relation between geometric albedo and po¬ 
larization properties for the asteroids is still an important 
task even in the era of large thermal radiometry surveys, 
whose results are valid in terms of statistical distribution 
among large samples of the population, but can be strongly 
inaccurate for what concerns single objects. 

It is therefore very important to optimize the perfor¬ 
mances of the polarimetric technique, as an effective tool to 
estimate the albedo of the objects, with particularly impor¬ 
tant applications to the physical characterization of newly- 
discovered, and potentially hazardous near-Earth objects. 
We stress again that, once a reliable calibration of the re¬ 
lation between albedo and polarimetric parameters is avail¬ 
able, the albedo values obtained by polarimetric data are 
not affected at all by uncertainties due to poor knowledge of 
the absolute magnitude of the objects, a relevant advantage 
over other possible techniques. 

In this paper we have carried out an extensive analysis, 
based on the idea of using for calibration purposes a still 
limited number of asteroids for which we can be reasonably 
confident to know the albedo with good accuracy. We have 
used this sample to obtain a new calibration of the classical 
slope - albedo law, with the polarimetric slope being derived 
from available data using different possible approaches. We 
have also analyzed other possible relations, including a new 
calibration of the classical Pmin - albe do relation, the mor e 
recently proposed p* - albedo relation dMasiero et al.l[2012tl , 
and a new relation based on the polarimetric parameter, 
introduced for the first time in this paper. The resulting 
values of the various polarimetric parameters for 22 asteroids 
considered in our analysis, and the corresponding values of 
albedo are given in Tables 2 and 3, respectively. 

The extensive analysis presented in this paper produced 
a variety of interesting results. For what concerns the “clas¬ 
sical” slope - albedo and Pmin - albedo relations, we conhrm 
that it is difficult to hnd a calibration which can ht accu¬ 
rately objects of all albedo classes. The results improve very 
much when the asteroids having low albedo, below 0.08, 
are not considered in the analysis. The presence of aster¬ 
oids having significantly different albedo but nearly identi¬ 
cal polarimetric slopes can even suggest that a parabolic ht 
is more suited to represent the h - albedo relation, as shown 
in Section Hr] (see Fig. IS]). 

The calibration of the h albedo and Pmin - albedo 
relations that we obtained by excluding low-albedo aster¬ 
oids from our analysis gives good results, with uncertain¬ 
ties within ±20% for medium- and high-albedo objects. For 
asteroids exhibiting a steep polarimetric slope, or having 
a polarimetric slope computed on the basis of only a few 
observations, the classical h - albedo relation can still be 
used, using the calibration coefficients found by considering 
all the available calibration objects, since the expected er¬ 
rors should be still acceptable. This is encouraging, because 
whenever there are at least a few measurements obtained 
at phase angles larger than 14°, it is generally possible to 


derive a value of the polarimetric slope h, and the resulting 
albedo values turn out to be reasonably accurate in abso¬ 
lute terms, although the relative error can be above 30% 
for low-albedo objects. As opposite, the use of Pmin should 
always be discouraged for asteroids having a deep negative 
polarization branch, as seen in Section FS.II 

If one wants to avoid the complication of using differ¬ 
ent calibrations of the slope - albedo relation for asteroids 
belonging to different albedo classes, other polarimetric pa¬ 
rameters can also be adopted. In particular, both and p* 
require to have at disposal phase - polarization curves of a 
sufficiently good quality to be ht by means of an exponential 
- linear representation (Eqjl)). This cannot be done when the 
number of polarimetric measurements is too small, or the 
observations are concentrated over a too limited interval of 
phase angles. While the uncertainty of albedo determina¬ 
tions based on the p* parameter seems to be reasonable but 
not really negligible, the most accurate albedo determina¬ 
tions can be obtained when the parameter can be reliably 
determined from the available data. The advantage of using 
4^ is that this parameter seems to be suited to give accurate 
values of albedo for both bright and dark asteroids. 

Of course, there are some caveats to be taken into ac¬ 
count. For instance, there are asteroids, like the so-called 
Barbarians, which exhibit pe culiar phase - polarization 
curves Isee ICellino et aP 120061 1. The determination of the 
geometric albedo for Barbarian asteroids is a problem, be¬ 
cause any derivation of the albedo using relations valid for 
the rest of the population may be misleading. However, it 
is also true that more data are needed to better understand 
the situation. 

In fact, the major problem of asteroid polarimetry today 
is certainly a serious lack of data. Further progress in this 
field will require the use of dedicated telescopes, to fill the 
gap with the amount of information already available from 
other observing techniques. 

Polarimetry has been so far severely under-appreciated 
as an essential tool for physical characterization of asteroids 
and also other classes of small bodies in the Solar System. 
The present paper summarizes the current state of the art 
for what concerns the calibration of the relation between 
polarimetric properties and albedo, and makes some further 
steps forward, with the introduction of the 4/ parameter, 
which seems to be a new useful tool to obtain reliable values 
of asteroid albedos. 

In a separate paper, to be submitted soon for publica¬ 
tion, we exploit the results of the present analysis to derive 
the al bedos of a fairly large number of objects for which we 
lack a IShevchenko and Tedescj ll2006l l determination, and 
we analyze the distributions of other parameters that char¬ 
acterize the phase - polarization curves of main belt aster¬ 
oids. 
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